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It has been shown experimentally that the ®Li nucleus shrinks by adding a A particle while 
its spectrum is also compressed. We discuss the compatibility of these two effects, contrasting 
also with a relation between the shrinkage effect and the stability of the spectrum in the ^Be 
hypernucleus. To this end, we employ two-body d-\}ie and a-^Ife cluster models for the \Li 
and aBc hypernuclei, respectively. We first argue that a Gaussian-like interaction between two 
clusters leads to a stabilization of the spectrum against an addition of a A particle, even though 
the intercluster distance is reduced. In the case of ^Li, the spin-orbit interaction between the 
intercluster motion and the deuteron spin has to be considered also. We show that the shrinkage 
effect makes the expectation value of the spin-orbit potential larger, lowering the excitation energy 
for the ''Li(3+) ® A(l/2+) level. 

PACS numbers: 21.80.+a,23.20.Lv,21.60.Gx,27.20.-f n 



I. INTRODUCTION 

One of the main interests in physics of hypernuclei is 
to investigate the impurity effect of A particle on the 
structure of atomic nuclei. This includes the change of 

deforma- 
m m, and 



e.g., nuclear sizeUj , the density distribution 
tion propertiesJSHlOj , the neutron drip- line 
fission barrier [13j. A characteristic feature of A particle 
is that it is free from the Pauli principle from nucleons, 
and thus it can deeply penetrate into the nuclear interior. 
It has been predicted that the A particle in the center 
of a nucleus attracts surrounding nucleons, leading to a 
shrinkage of nuclear size [l| . 

High precision 7-ray spectroscopy measurements have 
by now been systematically carried out for the p-shell A 
hypernuclei [14]. Such experiments have revealed that 
the electronic quadrupole transition probability, B{E2), 
from the first excited state (3+) to the ground state (1+) 
of ^Li is considerably reduced when a A particle is added 
p^ . By adopting a simple a-d and ^He-d cluster mod- 
els for ^Li and ^Li, respectively, the reduction in B{E2) 
has been interpreted as a reduction of the intercluster 
distance by 19%[lB|, in accordance with theoretical cal- 
culations [ll.[l6j. 

Although this interpretation has been well accepted, 
the spectrum of ^Li shows a somewhat puzzling feature. 
That is, one would naively think that the reduction of 
the intercluster size leads to a smaller moment of in- 
ertia, and thus the rotational excitation energy would 
increase. This has indeed been observed in the ^^C 
hypernucleus [l3| (notice that only the energy of the 3/2+ 
state has been measured in Ref. [l3|- However, the the- 
oretical calculation suggests that the partner state of the 
doublet (5/2+) lies below the 3/2+ state only by at most 
0.36 MeV [III- This corresponds to the spin averaged 
energy of the doublet of 4.66 MeV, that is compared 
to the energy of the 2+ state in ^^C, 4.44 MeV). On 
the contrary, the observed difference in the spin aver- 
aged energy between the first (l/2+,3/2+) and the sec- 



ond (5/2+, 7/2+) doublet levels in ^Li is 1.858 MeV 
that is smaller than the excitation energy (2.186 MeV) 
of the first excited state of ^Li (see Fig. 1). For the ^Be 
hypernucleus, the experimental data show that the spin 
averaged energy of the *Be(2+)(g)A(l/2+) doublet states 
(3/2+, 5/2+) is almost the same as the energy of the first 
2+ state of the *Be nucleus [l^, [2l[, despite a similar 
shrinkage effect, expected as in the ^Li nucleus[l|. The 
stability of the rotational energy of ^Be was shown to re- 
main the same even when the A particle is replaced by a 
charmed baryon Ac [l^l. Apparently, these behaviors of 
the spectrum of the ^Li and ^Be cannot be understood 
in analogy of e.g., a classical rotor of diatomic molecules. 
We summarize the experimental spaectra in Fig. 1. 



In this paper, we clarify the relation between the two 
contradictory effects of A particle on the structure of light 
atomic nuclei, that is, the shrinkage of the size and the 
compression/stabilization of the spectrum. An impor- 
tant point to remember is that the a cluster structure is 
well developed in the ground state of ^Li and ^Be, while 
the ground state takes the mean-field-like structure in 
^^C and in heavier nuclei. Our aim in this paper is thus 
to discuss how the spectra of ^Li and ^Be are compati- 
ble with the shrinkage of their size from a viewpoint of a 
two-body cluster model. 



The paper is organized as follows. In Sec. II, we first 
discuss the impurity effect of A in the ^Be nucleus. Using 
the a-^He cluster model together with the semi-classical 
Bohr-Sommerfeld quantization rule, we investigate how 
a A particle influences the rotational spectrum. In Sec. 
Ill, we discuss the rotational spectrum of the ^Li nu- 
cleus. We show that the spin-orbit interaction due to 
the deuteron spin plays an important role in lowering 
the excitation energy of the ^Li(g)A(l/2+) level. We then 
summarize the paper in Sec. IV. 
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FIG. 1: The experimental low-lying spectra for ^Li, "aLi, 
**Be, iBe, ^^C, and ^iC nuclei. The energies are denoted 
in the units of MeV. The energies for the l"'"(g)l/2^ and 
2+®l/2"'" levels of the ^Li nucleus are obtained by spin- 
averaging the observed level energies for the 1/2+ (0 MeV) 
and the 3/2+ (0.069 MeV) states, and the 5/2+ (2.05 MeV) 
and the 7/2+ (2.52 MeV) states 19], respectively. The energy 
for the 2+®l/2+ level of the ^Be is the spin averaged energy 
between the 5/2+ (3.024 MeV) and 3/2+ (3.067 MeV) states 
[2H . In the nucleus, only the energy for the 3/2+ state 
has been measured [l3|, while the energy for the 5/2+ state 
is expected to be lower than the energy of the 3/2+ state by 
at most 0.36 MeV \M- 



II. STABILIZATION OF LEVEL SCHEME FOR 

aBE 

We first discuss the level scheme of the ^Be hypernu- 
cleus, and clariiy the reasoning why the influence of the A 
particle in the level spacing is negligibly small. For this 
purpose, we use a simple two-body cluster model with 
a-^He configuration. Assuming a local potential, an in- 
tercluster potential at an intercluster distance R reads 

V^^MR) = VlaHR) + V^^\R) + KA(i?), (1) 



where V^^'' 



and Vaa^ 



are the nuclear and the Coulomb 
parts of the potential between the a particle and the core 
nucleus (^He) of the ^He, and VaA is the potential be- 
tween the A particle in ^He and the a particle. Assuming 
that the center of mass of ^He is identical to the center of 
mass of the core nucleus ^He, we use the same potential 



AN) 
^ aa 



and Vaa^ . That is, 



as that given in Ref. [23[ for Va 

with Fo=122.6225 MeV and a=0.22 fm~^ and 
l^i^ni?) = ^erf(/3i?2), 



(2) 



(3) 



with /3=0.75 fm ^, where erf(x) is the error function. For 
Vqa , we assume that it is given by the direct part of the 



double folding potential, that is, 

VaK{R) = J drNdrAPa{rN)pA{rA)vNA{R + rN~rf,), 

(4) 

where vna is a nucleon-A particle interaction, while pa 
and Pa are the density distributions for the a particle and 
the A particle in ^He, respectively. We assume Gaussian 
density distributions for pa and pA, 

p^{r) = 4(7r6„)-3/2e-'-Vb^„^ (5) 
PA(r) = (7r6A)-3/2e--^/6i. (6) 

Following Ref. flj, we take ba — 1.358 fm, and 6a = 
^(4MAr + A'/a)/4Ma 6a, where Mn and Ma are the 
mass of nucleon and A particle, respectively. For ^ata, 
we take the central part of the potential given in Ref. 



VNAir) = vo e 



(7) 



with 6,, 



1.034 fm, but we adjust the strength vq so 

9- 



as to reproduce the energy of the ground state of ^Be 
from the threshold of the a-l-^He configuration (that is, 
—3.50 MeV). Since pa, Pa, and vna are all given in a 
Gaussian form, the double folding integral in Eq. (|3]) can 
be evaluated analytically as 



bl + bl + bl 



3/2 



exp 



bl + bl + hi 



(8) 

With the intercluster potential so constructed, we solve 
the Schrodinger equation for the relative motion for each 
angular momentum L, 



sUi(i?)=0, (9) 



mA)/(8Mjv + Ma) is the 
reduced mass for the relative motion between a and a He. 
The Pauli principle is taken into account by excluding 
those states which satisfy 2n + L < A i23| , where n is the 
radial node of the wave function. 

The E2 transition probability can be computed with 
the intercluster wave functions iPl{R). In the two-body 
cluster model of {Ai,Zi) + (^2,^2), the E2 transition 
operator reads 




where 



eE2 



Te2 = eE2i?'r2^(i?), 



Ml + M2 



Ml + M2 



(10) 



Zi, (11) 



is the E2 effective charge. Mi and M2 being the mass of 
the fragments 1 and 2, respectively. 

The calculated spectra for the ^Be and ^Be are 
shown in Fig. 2. Since we do not include the spin- 
dependent part of nucleon-A interaction, the doublet 
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FIG. 2: Low-lying spectra for *Be and aBb obtained with 
the two-body cluster model of a+a and a+%lle, respectively. 
The energies are denoted in the units of MeV. The energies 
measured from the two-body thresholds are also shown in the 
parentheses. The arrows indicate the reduced E2 transition 
probabilities, B{E2), from the first excited state to the ground 
state. 



states (3/2+,5/2+) of «Be(2+) A(l/2+) are degener- 
ate in energy. In order to obtain the energy and the 
wave function for resonance states, we use a bound state 
approximation. That is, we have replaced the interclus- 
ter potential for R larger than the barrier radius with a 
constant so that the wave functions are confined inside 
the potential barrier [24- 2()j. We have confirmed that the 
energies obtained with this procedure are close to those 
energies that give the maximum of the energy derivative 
of the scattering phase shift. For instance, for the 2+ 
state of ^Be, the energy derivative of the phase shift is 
maximum at E=2.80 MeV, while the energy obtained in 
the bound state approximation is 2.85 MeV. 

From Fig. 2, one sees that the excitation energy 
slightly increases due to the addition of A particle. We 
have confirmed numerically that the potential which we 
use well reproduces the experimental phase shift for a+a 
scattering |23|]. If the resonance states were properly anal- 
ysed, this potential would therefore yield the excitation 
energy of the 2+ state in ®Be of 3.04 MeV (see Fig. 1). 
Then, the energy shift due to the addition of A particle 
is estimated to be small indeed, that is, as small as 0.05 
MeV. In contrast to the excitation energies, the B(E2) 
value from the first excited state to the ground state is 
altered drastically. The large change of B{E2) value is 
due to the fact that the attraction of the A particle makes 
the resonance levels of ^Be turn to the bound states in 
^Be. As a consequence, the wave functions become spa- 
tially much more compact, as shown in Fig. 3. The 
root-mean-square distance between the fragments is 4.90 
fm and 3.24 fm for the L — Q states of ®Be and ^Be, 
respectively. The energy change due to the A particle 
shown in Fig. 2, that is, 0.14 MeV, is much smaller than 
what would have been expected from a classical rotor. 
El — L{L + l)h^ /2iiR^ , which corresponds to the en- 
ergy shift of 2.64 MeV. Notice that the B(E2) value for 
the transition from the 2+ to the 0+ states in *Be is 
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FIG. 3: (Color online) The radial wave functions ul{R) de- 
fined as ^l{R) = ul{R)/R ■ Ylm{R) obtained with the two- 
body cluster model. The dashed and the solid lines are for 
the *Be and aBb nuclei, respectively. The upper panel shows 
the wave functions for the ground state with L = 0, while 
the lower panel shows the wave functions for the first excited 
state with L — 2. 



consistent with the value obtained in Ref. [27| using a 
damping factor method to evaluate the integral with the 
resonance wave functions, although this value is^ rather 
large as compared with that obtained in Ref. [Ij using 
the harmonic oscillator expansion of the wave functions. 

We now ask a question why the spectrum of ^Be is 
not influenced much, despite the fact that the wave func- 
tions are considerably altered. In order to address this 
question, the semi-classical approximation may be useful. 
In Ref. [23|, Rowley used the Bohr-Sommerfeld quanti- 
zation rule to show that the spectrum is approximately 
given by 



El — El=o 



^X(L- 



(12) 



for a Gaussian intercluster potential, V{R) — Vq e~^^ . 
We numerically confirm this relation in Appendix A. This 
formula indicates that the spectrum depends on neither 
the depth parameter Vq nor the intercluster distance, but 
only on the range parameter /3 of the potential as well as 
the reduced mass /i. 

We show the intercluster potentials for the present 
problem in Fig. 4 by the solid lines. The potential for the 
a-|-^He is deeper than that for a + a by 15.6 MeV, due to 
an additional attraction caused by the A particle. These 
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FIG. 4: (Color online) The intercluster potential for the a + a 
system (the upper panel) and for the Q+^He system (the 
lower panel). The dashed lines show the result of the fitting 
with a single Gaussian function. 



potentials are actually well fitted with the Gaussian form 
with Vo = -117.8 MeV and /3 = 0.222 fm-^ for ^Be and 
with Vo = -131.2 MeV and /3 = 0.250 fm^^ for ^Be. 
The quality of the fit can be seen in Fig. 4. The reduced 
masses, on the other hand, are 1877.8 MeV and 2120.6 
MeV for the a + a and the a+^He systems, respectively. 
Here, we have used the average value of the proton and 
neutron masses as the nucleon mass. The factor jSh^/Sn 
in Eq. dH]) is thus 0.574 MeV for the a + a system and 
0.572 MeV for the the a+^He system. It is remarkable 
that these values are considerably close to each other. 
Evidently, a Gaussian-like intercluster potential has a re- 
sponsibility to make the spectra resemble to each other 
between *Be and a Be, despite that the absolute value of 
the energies and thus the radial dependence of the wave 
functions are considerably changed due to an addition of 
a A particle. 



III. LEVEL SCHEME FOR ^LI: ROLE OF 
SPIN-ORBIT INTERACTION 

Let us discuss next the spectrum of the ^Li nucleus. If 
we assume that the nucleus ^Li takes the a-l-deuteron (d) 
structure, a big difference of this nucleus from *Be=a-|-a 
is that one of the fragments {i.e., the deuteron) has a 
finite spin {S = 1). The potential between a and d then 



reads [29|-l31|, 



Va^diR) = VoiR) + Vi{R)L-S + V2{R) 



i?2 



(13) 

The spectrum of the ground rotational band of ^Li can 
be understood as follows j31i]. This band arises from the 
states with 2n -|- L = 2, that is, L = (Is) and L — 2 
(Od). By combining with the deuteron spin S' = 1, the 
Is state yields the total spin and parity of /'^=1+, while 
the Od state yields /'^ = 1"'",2+, and 3+ triplet states. 
These triplet states are degenerate in energy in the cen- 
tral potential Vo{R), but they split by the spin-orbit in- 
teraction Vi{R)L ■ S as well as by the tensor interaction 

V2{R) [(5 • - 5'V3 . 

In order to investigate how the spin-orbit interclus- 
ter potential affects the spectrum of the ^Li nucleus, we 
closely follow Ref. [3l| and use the semi-microscopic clus- 
ter model of Buck et aZ. js^]- To do so, we assume that the 
^Li nucleus takes the d-|-^He structure. In this model, 
the intercluster potentials Vo{R),Vi{R), and V2{R) are 
constructed based on a core-|-p-|-n three-body model. We 
regard the a particle and ^He as the core nucleus for ^Li 
and ^Li, respectively. In the following, for simplicity, 
we assume that the deuteron is in a pure s-state, and 
thus the tensor part of the intercluster potential van- 
ishes. This term plays an essential role in reproducing 
the quadrupole moment of ^Li^29-.31j. but its infiuence is 
expected to be small for the spectrum. With this approx- 
imation, the central and the spin-orbit potentials read 



VoiR) = drUo{r,R)x{rY 
Vi{R) = 



(14) 



2M, 



2{Mn + M,) 
dr 



Wo{r,R) + —Wi{r,R)\ xW^,15) 



where Mc is the mass of the core nucleus {i.e., "^He or 
^He) and x(r) is a s-state wave function for relative mo- 
tion between p and n in the deuteron. Uq, Wq and Wi are 
multipole components of the nucleon-core interaction. 



cN 



'U{R±r/2) - W{R±r/2)l- s, 



(16) 



where r is the relative distance between p and n in the 
deuteron, I is the relative angular momentum between 
the core and the nucleon, and s is the nucleon spin. That 
is, Uo, Wq and Wi are defined as 

U{R±r/2) = 47r^C/A(i?,r)n(±r)-n(-R),(17) 

A 

W{R±r/2) = An J2 Wx{R, r)Yx{±r) ■Yx{R). {18) 

X 

In the calculations shown below, following Ref. [sH , we 
employ the Hulthen form of the deuteron wave function. 
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that is, x(r) = •\/2ae~"'" with a=0.2316 fm"-'^, and the 
same potential as in Ref. (3T| (with a slight readjustment 
for the depth parameters) for the nucleon-a potential. 
The contribution of the A particle to the ^He-nucleon 
potential is estimated with a single folding potential, 



UNAir) = J drAPAirA)vNAir + rN -va), 



(19) 



where the A particle density pA and the nucleon-A inter- 
action vna are given in the previous section. As in Eq. 
(HI, this potential can be computed analytically as 



UNA{r) 



Wo 



bl + bl 



3/2 



exp 



We adjust the value of vo in order 
ground state energy of the ^ 



to 



bl 



(20) 



^Li nucleus (that is. 



MeV from the threshold of %lle+d) 
potential for ^He-nucleon is given by U (r 



reproduce the 
3.48 

The total central 
UNa{r) + 

UNAir), where UNa(r) is the central part of the nucleon- 
a potential. Since we neglect the spin-orbit interaction 
in the nucleon-A interaction, the A particle does not con- 
tribute to the spin-orbit part of the intercluster potential, 
Vi{R), except for the trivial mass factor in Eq. p3|) . 
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FIG. 6: (Color online) Low-lying spectra for Li and yvLi ob- 
tained with the two-body cluster model of Q-fd and ^He-l-d, 
respectively. The results obtained with and without the spin- 
obit intercluster potential are compared. The energies are 
denoted in the units of MeV. The energies measured from 
the two-body thresholds are also shown in the parentheses. 
The arrows indicate the reduced E2 transition probabilities, 
B{E2), from the first excited state to the ground state, where 
the calculated values are denoted in the units of e^im^. 
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FIG. 7: (Color online) (a) The radial wave functions for the 
ground state with L — for ^Li (the dashed line) and ^Li 
(the solid line), (b) Those for the first excited state with 
L = 2 obtained without including the spin-orbit intercluster 
potential, (c) The overlap between the spin-orbit potential 
and the unperturbed wave functions for L = 2. (d) The same 
as Fig. 7(b), but those obtained with including the spin-orbit 
potential. 



FIG. 5: (Color online) The central part (the upper panel) and 
the spin-orbit part (the lower panel) of the intercluster poten- 
tial. The dashed lines show the potentials for the ^Li nucleus, 
while the solid lines show those for the ^Li hypernucleus. 

The central part Vo{R) and the spin-orbit part Vi{R) 
of the intercluster potentials are shown in Fig. 5. For 
the central part, we also include the Coulomb potential 
evaluated in the same way as in Ref. 'sT] . The solid and 
the dashed lines are the potentials for the \lle-d and the 



a-d, respectively. The depth of the central part increases 
by 5.98 MeV due to an addition of a A particle, while the 
change of the spin-orbit potential is much smaller. The 
spectra obtained with these potentials are shown in Fig. 
6. In the absence of the spin-orbit potential, Vi{R), an 
addition of the A particle slightly increases the energy 
of the first excited state, from 3.72 MeV in ^Li to 3.80 
MeV in ^Li. This behavior is similar to the change of 
the spectrum of ^Be shown in Fig. 2. The radial wave 
functions for the ground state and the first excited state 
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obtained without the spin-orbit potential are shown in 
Figs. 7(a) and 7(b), respectively. As in the case of a Be, 
the shrinkage effect of the A particle is substantial [ll.[l6l]. 

Let us now consider the effect of the spin-orbit po- 
tential. We fist treat it by the first order perturbation 
theory, that is, the energy shift for L = 2 and / = 3 is 
approximately given by 



ASl=2,/=3-2 / dRViiR) u^Ll^iR) 



(21) 



Here, the factor 2 is the eigenvalue of the operator L ■ S 

for L = 2 and / = 3, and ^^^^2 ^^'^ unperturbed wave 
function for L — 2 obtained without the spin-orbit po- 
tential. Using the wave functions shown in Fig. 7(b), 
we obtain AEl=2.i=3 = -1-22 MeV for ^Li and -1.84 
MeV for ^Li. Fig. 7(c) show the integrand of Eq. (|2T]) . 
One can clearly see that the overlap of the wave function 
with the spin-orbit potential shown in the lower panel 
of Fig. 5 increases significantly for the ^Li hypernu- 
cleus due to the shrinkage effect, leading to the larger 
value of the energy shift AEl=2,i=3- The spectra ob- 
tained by treating the spin-orbit potential exactly, and 
the corresponding wave functions are shown in Figs. 6 
and 7(d), respectively. The effect of the spin-orbit po- 
tential is indeed much larger in ^Li due to the shrinkage 
effect, and the energy of the "^0(3+) <S) A(l/2+) state 
appears lower than that of the ^Li(3"*') state. We thus 
conclude that the spin-orbit potential due to the finite 
deuteron spin plays an important role in compressing the 
spectrum of the ^Li hypernucleus. We mention that the 
calculated B{E2) value from the 5/2+ member of the 
^Li(3+) (g) A(l/2+) state to the 1/2+ state in the ground 
state <5Li(l+) ® A(l/2+) doublet is 3.03 e^fm"*, that is in 
good agreement with the experimental value, 3.6±0.5tQ'4 
e^fin^ [l^. The calculation also reproduces well the ex- 
perimental value for the energy of the ^Li(3+) A(l/2+) 
state, 1.86 MeV (See Fig. 1). 



IV. SUMMARY 

We have investigated the structure of ^Be and ^Li hy- 
pernuclei using two-body cluster models of a + d and 
^He-l-fi, respectively. We particularly discussed the rela- 
tion between the shrinkage effect and the spectra of these 
hypernuclei. We first showed that the two-body cluster 
model yields only a small change in the spectra between 
^Be and ^Be even though the intercluster distance is sig- 
nificantly reduced due to an addition of a A particle. We 
argued based on the Bohr-Sommerfeld quantization rule 
that this is caused by the fact that intercluster potentials 
can be well fitted with a single Gaussian function, for 
which the spectra are approximately independent of the 
depth of the potential. A similar effect would be expected 
also for ^Li if the deuteron were a spin-less particle. In 
reality, a deuteron has a finite spin (S'=l), and the spin- 
orbit potential is present in the intercluster potential. 
We have shown that the shrinkage effect leads to a large 



overlap between the wave function for the excited state 
and the spin-orbit potential for the ^Li hypernucleus as 
compared with ^Li. As a consequence, the spin-orbit po- 
tential acts effectively larger in ^Li, eventually leading to 
the compression of the spectrum. 

These behaviors of the spectra are peculiarities of a 
two-body cluster structure. In this sense, ^Li and ^Be 
are exceptional cases, as the ground states exhibit a 
well-developed a cluster structure. For heavier nuclei, 
the cluster structure appears in excited states while the 
ground state takes a mean-field-like configuration. In 
such situations, the shrinkage of radius would push up 
the spectrum, as has been seen experimentally in "'^^C. 
The deformation degree of freedom also comes into a play 
there. An interesting case for a future investigation is 
the hypernucleus ^^F, in which the mean-field structure 
may be mixed with the ^J^O+d cluster structure. A 7- 
ray spectroscopy measurement has been planed on this 
hypernucleus at the J-PARC facility [33,, i3^] . This will 
be the first for the sd-shell hypernuclei. The ground and 
the first excited states of the core ^®F nucleus has a spin 
and parity of 1+ and 3+, respectively, which are the same 
for ^Li. Thus the measurement of the i*^F(3+)(g)A(l/2+) 
level will provide a strong stimulus to further studies on 
the sd-shell hypernuclei and beyond. 
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Appendix A: Rotational spectrum of a Gaussian 
potential 

In this Appendix, we numerically check the perfor- 
mance of the approximate formula given in Eq. (jl2[) 
for the rotational spectrum for a Gaussian potential, 
V{R) = Vq e~^^ . We consider here the a + a system, 
and use the same geometry for the potential as that given 
in Eq. i.e., jS — 0.22 fm~^. In this Appendix, we do 
not include the Coulomb interaction between the two a 
particles. Fig. 8(a) shows the energy difference between 
the L — 2 state and the L — state for 2n+L — 4. Notice 
that without the Coulomb interaction the L = 2 state is 
bound only for Vq < —131.5 MeV. The solid line shows 
the numerical result of the Schrodinger equation, while 
the dashed line is obtained with Eq. (fT2|) . One can see 
that Eq. works well, the deviation from the exact re- 
sult being about ±0.2 MeV for this system. The energy 
difference indeed depends on the depth parameter only 
weakly, while the energy and the root mean square (rms) 
distance between the two a particles for the L = state 
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varies from -3.269 MeV and 3.58 fm at Vq = -132 MeV 
to -69.538 MeV and 1.93 fm at Vq = -300 MeV, respec- 
tively (see Figs. 8(b) and 8(c)). As sliown by the dotted 
line in Fig. 8(a), this is in marked contrast with the 
energy difference for a classical rotor, AE — 6ft^/2/ii?^, 
where /i — Ma/ 2 is the reduced mass for the a + a system 
and R is the rms intercluster distance. 
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FIG. 8: (Color online) (a) The energy difference between the 
L = 2 (Id) and the L = (2s) states for a Gaussian potential 
V{R) — Vo e~^^ for a a + a system as a function of the depth 
parameter — Vq. The solid line is obtained by numerically 
solving the Schrodinger equation, while the dashed line is the 
result of the approximate formula given in Eq. (|12p based 
on the Bohr-Sommerfeld quantization rule. The dotted line 
denotes the energy difference for a classical rotor, evaluated 
with the root-mean-square (rms) intercluster distance for the 
2s state shown in Fig. 8(c). (b) The energy of the L = 
(2s) state as a function of the depth parameter — Vo for the 
Gaussian potential, (c) The rms intercluster distance for the 
2s state. 
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